Measuring the phonon-assisted spectral function by using a non-equilibrium 

three-terminal single-molecular device 
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The electron transport through a three-terminal single-molecular transistor (SMT) is theoretically 
studied. We find that the differential conductance of the third and weakly coupled terminal versus 
its voltage matches well with the spectral function versus the energy when certain conditions are 
met. Particularly, this excellent matching is maintained even for complicated structure of the 
phonon-assisted side peaks. Thus, this device offers an experimental approach to explore the shape 
of the phonon-assisted spectral function in detail. In addition we discuss the conditions of a perfect 
matching. The results show that at low temperatures the matching survives regardless of the bias 
and the energy levels of the SMT. However, at high temperatures, the matching is destroyed. 

PACS numbers: 73.63.Kv, 71.38.-k, 85.65.+h 
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I. INTRODUCTION 

In the past decade, transport properties of single- 
molecule transistors (SMT) have attracted great atten- 
tion due to the potential application for the new gen- 
eration of electron devices. Because of the intrinsic vi- 
brational freedom in molecules, the molecular electronic 
transistor also provides a new opportunity for exploring 
the vibration-electronic (i.e. electron-phonon) interac- 
tions at single molecule level. The electron-phonon inter- 
action in a SMT leads to some interesting effects, such 
as phonon-assisted tunneling, the red shift of SMT en- 
ergy levels, and generation of the thermal energy. Such 
features are interesting and have been extensively inves- 
tigated both experimentally and theoretically in recent 
years. Park et alX experimentally studied the current- 
bias (I-V) characteristic of an individual Cgo molecule 
connected to gold electrodes and have observed the ob- 
vious phonon-assisted tunneling sub-steps in the I-V 
curves. Later the current of a suspended individual 
single- wall nanotube device is measured and two phonon- 
assisted sub-peaks on the two sides of the main resonant 
peak are clearly visible in the differential conductance 
versus the gate voltage, which is due to the radial breath- 
ing phonon mode^i 3 . Very recently, also in the device of 
a suspended single-wall nanotube but with much lower 
temperatures, the higher order phonon-assisted sub-steps 
on the I-V curves have been experimentally demonstrated 
by Sapmaz et alA 

On the theoretical side, there have also been a large 
amount of studies on the quantum transport behavior of 
a SMT or a quantum dot (QD) coupled to local phonon 
modes ,£ i £iL£ i 9 i 10 i 11 About 15 years ago, Wingreen et al. 
studied the electron transport through a QD coupled to 
the phonon modes by combining the scattering theory 
and the Green's function method, and the phonon-inducc 
transmission sidebands were found.— Using the real-time 
renormalization-group method, Keil et ai£ have investi- 
gated the quantum transport phenomena through cou- 



pled QDs with a phonon bath, and a solution for station- 
ary current is obtained. In addition, the shot-noise spec- 
troscopy of the current of a SMT having a local phonon 
mode is reported by Zhu and Balatsky. 8 Since the cur- 
rent, the conductance, the shot-noise, etc., are all closely 
related to the local electronic spectral functions A(u>), the 
spectral function has also been extensively studied. The 
spectral function is found to be strongly dependent on 
the positions of the SMT (or QD) energy levels and the 
tunneling strengths T between the leads and the SMT 
(or QD). 9 On the other hand, by using a different ap- 
proximation, Chen et al. also investigated the spectral 
function and the current through a SMTji and the spec- 
tral function in their results exhibits intriguing features. 
At low temperatures and under the condition that the en- 
ergy level is near the Fermi surface, the side peaks in the 
spectral function are clearly non-Lorenzian in shape. The 
side peaks on one side changes gradually while the other 
side change abruptly with changing of energy. These are 
quite different comparing with the previous results. Due 
to the importance of the spectral function A(u>) as well as 
that the spectral function of a SMT coupled to a phonon 
mode is not well understood (i.e. the results of A{u) 
are qualitatively different for using different approxima- 
tions), thus it is quite beneficial to design an experimen- 
tal set-up to directly measure the spectral function of a 
SMT. 

Recently, some studies have investigated the Kondo ef- 
fect in QDs coupled to the three terminals, in which the 
third terminal acts as an exploring tip to measure the 
Kondo peaks in the spectral function] 12 : 13 Can the in- 
triguing characteristics of the phonon-assisted side peaks 
in the spectral function also be explored in a three- 
terminal set-up? It is the purpose of this work to theo- 
retically analyze the feasibility of this scheme. We con- 
sider the system of either SMT or QD coupled to three 
leads and a local phonon mode. Here the third lead is 
introduced as an exploring tip. We find that the spectral 
function quite often matches the differential conductance 
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of the third lead. If the third lead is weakly coupled and 
the temperature T is low {kgT <C Hu>q with ojq being the 
phonon frequency), this matching is almost perfect (in- 
cluding the abrupt changes associated with the phonon- 
assisted side peaks), so the spectral function, in particu- 
larly the phonon- assisted side peaks, can be directly de- 
tected using the differential conductance. On the other 
hand, if the coupling of the third lead is large but still 
at low temperature, the matching is in qualitative agree- 
ment. However, at high temperature, the matching is 
destroyed even for weakly-coupled third lead. 

The rest of this paper is organized as follows. We in- 
troduce the model in Sec. II and derive formulas of the 
spectral function and the differential conductance in Sec. 
III. In Sec. IV, we show the numerical results and present 
discussions of those results. Finally, a brief summary is 
given in Sec. V. 

II. MODEL HAMILTONIAN 

The device under our consideration is illustrated in 
Fig.l. It consists of a SMT or QD connected to three 
metallic leads. An electron in the SMT is also coupled 
to a single phonon mode. Due to the big energy gap 
between two quantum levels in the SMT, only one rel- 
evant quantum level is considered. The Hamiltonian of 
the device is written as: 

H = H Leads + Hph + Hn + Ht ■ (1) 

The first two terms are, respectively, the Hamiltonian for 
electrons in the three leads and the Hamiltonian for the 
phonon part, 

Hleads = ^ ^ £ ak^ ak ^ak ; (2) 
ct,k 

Hph = uJo^b. (3) 

Here a — L, R, and 3 respectively represent the left, 
the right, and the third leads, and c ak (c ak ) creates (an- 
nihilates) an electron with the energy e ak in the lead 
a. Analogously, ^(b) is the phonon creation (annihila- 
tion) operator and ujq is the vibrational frequency of the 
phonon. The third term in Eq.(l) is 

H D = [e + \(b f + b)]dU, (4) 

where (d) is the electron creation (annihilation) opera- 
tor in the SMT with the energy level £o, and A describes 
the coupling strength between the SMT and the local 
phonon mode. The last term in Eq.(l) descries the tun- 
nelling coupling between the SMT and the three leads, 

H T = J2[V<* k cl k d + H.c}. (5) 

a,k 

It is often useful to take a canonical transformation 
with^^i? = e s He~ s and s = {\/u ){tf -b)dU. Under 



this canonical transformation the Hamiltonian (1) varies 
into 

H = H e i + Hph, (6) 

where 

Ha ^^e ak c^ ak c ak +e dU + ^2{V ak c ak d + H.c\, (7) 
and 

H ph = Luotfb. (8) 

Obviously, due to the electron-phonon interaction, the 
energy level £o of the SMT is renormalized to Eq — 
Eq — gu)Q with g = (A/wo) 2 , an d the tunnelling ma- 
trix element V k is varied into V k = V k X where X = 
exp [— (A/cl>o)(&' + b)]. Up till now, no approximation has 
been made, the Hamiltonian (6) is completely equivalent 
to the Hamiltonian (1)^ 

III. THE SPECTRAL FUNCTION AND THE 
DIFFERENTIAL CONDUCTANCE 

In this section, we calculate the spectral function A(u>) 
of the SMT and the differential conductance of the third 
lead. From the results by Meir and Wingreen, the spec- 
tral function A(u) and the current can be represented by 
the Green functions of the SMT as>±^ 

A{uj) — i{G > (uj) — G < (lj)] — i[G r {uj) — G a (uj)]. (9) 

Ja = \ J ^{r^M + ^H/aH]}. (io) 

Here f a (uj) = l/{exp[(uj — /i Q )/fcsT] + 1} is the Fermi 
distribution function with the chemical potential /j, a , 
r a (w) = 2ir^ k \V ak \ 2 8{e ak — u) describes the cou- 
pling strength between the lead a and the SMT, and 
Q<,>,r,a are standard lesser, greater, retarded, ad- 
vanced Green functions j 1 ^ '— Because of the existence of 
the electron-phonon interaction, it is difficult to directly 
solve these Green functions from the equation of mo- 
tion technique or the Dyson equations. As it is done 
in some previous papers?^ here we take the same ap- 
proximation to replace the operator X and X' in Hamil- 
tonian (6) with their expectation value (X) = (X') = 
exp [—g(N p h + 1/2)], where N p h = l/[exp(/3w )-l] is the 
phonon population^ This approximation is valid when 
the tunneling strengths r a are smaller than the electron- 
phonon interaction, i.e, T a -C A. After this approxi- 
mation, the Hamiltonian (6) is decoupled into two inde- 
pendent parts, electronic part and phonon part. Next, 
we also need to decouple the Green functions. In many 
previous papers^ they decouple the retarded (advanced) 
Green functions G r ' a directly. However, such decoupling 
has some defects as pointed out in a recent work by Chen 
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et aiJ£ Here we employ the decoupling method in the 
Ref.(— ), to directly decouple the lesser and greater Green 
functions G <,> instead of the retarded and advanced 
Green functions G r,a . After the decoupling, the lesser 
and greater Green functions G <> are: 



G<{t,t') = G^M') e-*(*'-*), 



(11) 
(12) 



where G (t, t') are the Green functions of the Hamil- 
tonian H e i, and $(i) is 

= g [N ph (l - e- ot ) + (N ph + 1)(1 - e -*"°*)] . (13) 

Using the identity^ e 2COsS = J2n=-Z I n (z)e mB , the 
greater and lesser Green functions can be expanded as 



G < (u>)= B n G (w + nua), 

n— — oo 
+ oo 



(14) 
(15) 



n— — oo 



where the coefficients B n = 

e -g(2N ph + l) e nu /2k B T In (2 gy /]Sr ph (N ph + 1) ) and 

I n (z) is the nth Bessel function of complex argument. 
Thus we can rewrite the spectral function A(ui) as 

A(u) = ^[G > (w) - G< (w)] 

+00 

= ^3 iB n [G (uj — nujo) — G (ui + nojo)I-16) 

n=— 00 

Following the standard derivation ) 14 i 16 ' 17 the self- 
energies £ of the coupling to the leads for the Hamil- 
tonian H e i can be easily obtained as: 

r {a \co) = ^|F Qfc | 2 5 l a) H = ^[A Q (c) T ^r tt (^]7) 

Q.fe a 

s < ( w ) = 53 |^| 2 <?<» = i£r Q H/ Q M, (is) 

Q.fe Q 

s>H = 53|y Qfc | 2 5 > fc H--z53r Q H[i-/ Q pa 



where r Q = r Q exp [— g(2N p h + 1)] since the tunneling 
elements V a k have been amended by electron-phonon in- 
teraction. To take the wideband limit ) 16 ' 17 i.e. to assume 
that r a then r a are independent with the energy w, the 
above self-energies reduce into: 



2 H 



(20) 
(21) 

(22) 



By using these self-energies, the dressed retarded (ad- 
vanced) Green function Cf" "\ then the dressed lesser and 

greater Green functions G ' can be readily obtained 
from Dyson equations and Keldysh equations: 

G r{a) (to) = ^(wj-f'"'^]- 1 , (23) 
G^w) = G r (u)E < ((j)G a (cj) = if{uj)A(u), (24) 
(^{u) = G r (a;)S > (a;)G a (a;) = -i[l - J{u)]A(d$>) 

where 



(a ) -e ) 2 + (rL + r fl , + r 3 ) 2 /4' 



(26) 



and 



7 , 1 = r L / L H + r fl / fl (a;) + r 3 /3M 
r L + r fl + r 3 

After solving G < ' > , the lesser and greater Green func- 
tions G <,> (w), the electronic spectral function A(uS) of 
the SMT, and then the current can be calculated from the 
above Eqs. (9), (10), (14), and (15), straightforwardly. 

At last, the differential conductance G 3 of the third 
terminal can be acquired by performing G 3 = 9J 3 /9V 3 

<*- E ^«./|{A(«)H-««)l 



n— — oo 



xp( W2 )/M+I( Wl )[i-7( Wl )]] 

-c[l - / 3 (w)]/sM[l - /3^3)P(wa) 
-c/aMAWtl - / 3 (^iP( Wl )}, (28) 

where cji = uj — nujQ, u>2 = u + nwo, and c = r 3 /(rz, + 

r fl + r 3 ). 



IV. NUMERICAL RESULTS AND 
DISCUSSIONS 

In this section, we study numerically the spectral func- 
tion A(cu) and the differential conductance G 3 . In the 
numerical calculation, the coupling strengths T L / R be- 
tween the left/right lead and the SMT is set to be unity 
(r^ = Tr = T = 1), as an energy unit. The main pur- 
pose in the present work is to study whether the curve 
of the spectral function A(u) versus the energy u> can 
map into the curve G 3 -V 3 , i.e. whether the intriguing 
phonon-assisted side peaks in the spectral function A{uj) 
can be observed by measuring the conductance G 3 of the 
third lead. First, let us show the spectral function of the 
two-terminal SMT with T 3 = 0. Note that this spectral 
function A{lo) is the object of our study. When the SMT 
is coupled to the phonon mode, one main characteristic of 
A(ui) is the appearance of the phonon-assisted side peaks. 
At the zero bias case {^l = I^r = 0) and the renormalized 
level Eq — 0, the side peaks are non-Lorenzian in shape, in 
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which one side of the side peaks still looks like the Loren- 
zian form but the other side drops abruptly (see Fig. 2). 
When f?o (i- e l^ol/r 3> 0) is far away from the chemical 
potentials the side peaks are asymmetry on the 

two sides of the main peak, and the side peaks disap- 
pear on one side. Furthermore, with a non-zero bias V 
(V = /il — Hr) or raising temperature T, these phonon- 
assisted side peaks exhibit more complex profiles. These 
characteristics of the spectral function A(u>) have been 
found in a previous study^ Our goal here is to propose 
a scheme to measure the spectral function A{ui) by using 
an extra third lead. 

We first study the zero bias case (/il = Mi?. = 0) with 
a weakly-coupled third lead = 0.01). Fig. 2 shows 
the differential conductance G3 versus the voltage V3 of 
the third lead for different renormalizcd level £q. For 
comparison, the spectral function A{lo) versus the energy 
oj for the two-terminal SMT device with r3 = is also 
shown in Fig. 2. When £0 = 0. Besides the main resonant 
tunneling peak, some extra phonon-assisted side peaks 
emerge in the curve G3-V3. The main peak is Lorenzian, 
but the side peaks exhibit non-Lorenzian characteristic. 
On the one side of the side peaks, the conductance G3 
falls abruptly from top to valley. In particular, we find 
that the curve of the conductance G3 versus V3 is in an 
excellent agreement with the curve of the spectral func- 
tion A(u>) versus the energy to. Not only are their side 
peaks located at the same positions, but also they have 
the same profiles. Even the abrupt drops overlap per- 
fectly. Thus, in this case by measuring the differential 
conductance G3, one obtains all information on the spec- 
tral function A(u>). Increasing £"0 to 1, the side peaks in 
the curve of G3-V3 are distributed asymmetrically on two 
sides of the main resonant peak. The right side peak is 
higher than the corresponding left side peak as shown in 
Fig. 2. However the curve of G3-V3 is still in an excellent 
matching with the spectral functions A(ui) versus w, even 
the complex structure of the first right side peak (about 
at V3 = 6) matches quite well. When £q is furthermore 
enhanced from 1 to 3, all side peaks are on the right of 
the main peak and all peaks are Lorenzian. Similarly the 
excellent agreement between the differential conductance 
and the spectral function arc still maintained. Combining 
above results we find that at the zero bias case, the curve 
G3-V3 is in excellent agreement with the curve A(uj)-lu re- 
gardless of the value of the level £0 (i-e- £o)- The reason 
lies in the fact that for the weakly-coupled lead, the trans- 
mission probability of the incoming electron with energy 
lo is mainly determined by the local spectral function 
A{u>). 

Next, we study the case with a finite bias V (V = 
— ^b,)- Fig. 3 shows the differential conductance G3 
and the spectral function A(lu) for /i£ = — /i_R = 3. In 
the finite bias V = 6 and the renormalized level £0 = 0, 
the phonon-assisted side peaks in the conductance G3 are 
symmetrically distributed on two sides of the main reso- 
nant peak, and the form of the side peaks are Lorenzian 
which is in contrast to zero bias case with non-Lorenzian 



side peaks (see Fig. 2). However, the characteristics of 
the spectral function are still reflected perfectly by the 
differential conductance, although the spectral functions 
for the zero bias and the non-zero bias have a large differ- 
ence. When £0 is reduced to —3, the left phonon-assisted 
side peak is obviously higher than the corresponding right 
side peak, the first left side peak becomes sharp, and the 
non-Lorenzian characteristic emerges again in the first 
right side peak. Although the spectral function A(uj) is 
so complex now, the curve of G3-V3 closely follows that 
of the spectral function, including the detail structure. 
With decreasing £ further, e.g. £0 = —6, the conduc- 
tance G3 has a large change. In the present case, the side 
peaks of the right hand disappear completely and the side 
peaks only exit on the left hand of the main peak. We 
find that the curve G3-V3 completely matches the curve 
A(uf)-u). In fact, for any bias V and any level £0 (i.e. £0), 
the curves G3-V3 and A{uj)-uo overlap perfectly regardless 
of complexity of the curve A(lj)-lu. This means that by 
measuring the differential conductance G3 of the third 
lead, the spectral function A(u>), including the intriguing 
characteristics due to coupled to the phonon mode, can 
be directly observed. 

In the above numerical investigation, the coupling T3 
between the SMT and the third lead is set to be rather 
weak (F3 = 0.01), and the temperature T is kept at very 
low (ksT = 0.02). For a strongly-coupled third lead 
and at high temperature, will the excellent agreement 
between the differential conductance G3 and the spectral 
function A(u) still survive? In this paragraph, we study 
the effect of the coupling strength T 3 . The temperature 
effect is addressed in the next paragraph. Fig. 4(a) and 
(b) show the conductance G3 and the spectral function 
A(uS) for different coupling strength T3. With increas- 
ing T3, the level in the SMT is widen, and the conduc- 
tance G3/G0 [G = (2e 2 /fr)(r 3 /r)] is overall reduced. 
To make a better comparison of the two curves G3-V3 
and A(u)-u>, we replace the conductance unit Go by an 
integral weighting factor G' and Fig. 4 (c) and (d) show 
G3/G0 versus V3, here G' is determined by the equa- 
tion f = f %±£kal. The results are as following: 
(i) When the third lead is weakly-coupled with a small 
^3 {F3 — 0.01), the differential conductance G3 and the 
spectral function A(u>) map into each other perfectly, as 
discussed before, (ii) With increasing coupling strength 
T3, the differential conductance deviates from the spec- 
tral function gradually. When T3 is in the same order of 
T (e.g T 3 = 0.4 or 1), the peak of the conductance G 3 be- 
comes lower and wider than that of the spectral function 
(see Fig. 4). But overall the conductance G3 still shows 
the similar profile of the spectral function, including the 
abrupt drop on the phonon-assisted side peak. In other 
words, the curve G3-V3 is still qualitatively the same as 
the curve A(u))-u> when T3 ~ T. Therefore, the probing 
lead is not necessarily need to be weakly coupled to the 
SMT, and the device can still work at T 3 ~ T. (iii) When 
T3 is much larger than T (e.g. T3 = 4r or more), the side 
peaks in the conductance G3 fade away, and the curves 
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of G3-V3 and A{uS)-u> have a large discrepancy. Since the 
properties [including the spectral function A(u)] of SMT 
are remarkably influenced when the coupling between the 
third lead (i.e. exploring terminal) and the SMT is very 
strong, the differential conductance can no longer exhibit 
the characteristics of the real spectral function. From 
Fig. 4, it can be concluded that the qualitative agreement 
between the curves of G3-V3 and A(lu)-ui is destroyed 

when r 3 > 2(r L + r R ). 

Now let us consider how temperature influence the re- 
lationship between the spectral function and the differen- 
tial conductance. As displayed in Fig. 5a, at low tempera- 
ture (T = 0.05r) the curve of the spectral function over- 
laps perfectly with that of the differential conductance, 
as being discussed before. When temperature is raised 
to O.fr, although the peaks of the differential conduc- 
tance are slightly lower (about a few percent) and wider 
than those of the spectral function, they still agree not 
only in their positions, but also in the shape of the peaks 
and the detail structure of the phonon-assisted side peaks 
(see Fig. 5b). With further increasing temperature T, the 
deviation between the conductance G3 and the spectral 
function A(lo) is more noticeable. When T = 0.5F, the 
peak heights of the conductance decrease to half of those 
of the spectral function. In particular, the shapes of the 
side peaks in the conductance and in the spectral function 
are clearly in disparity. The shapes of the side peaks of 
the spectral function are still asymmetric and the change 
is quite abrupt on one side of the side peak, but the 
shapes of the side peaks of the differential conductance 
are Lorenzian and symmetric (see Fig. 5c). While the 
temperature is equal to or larger than T, all side peaks of 
the differential conductance fade away, and the differen- 
tial conductance is no longer providing any information 
on the spectral function (see Fig.5d). Therefore, it is 
feasible only at the low temperature to observe the spec- 
tral function by measuring the differential conductance 
of the third terminal. In the Fig. 5, we choose zero bias 
and the renormalized level £0 located at zero. In fact, all 
conclusions remain for non-zero bias and any value of So ■ 

Unlike the coupling strength T3, the temperature 
strongly influences the comparability of the curves G3-V3 
and A(u))-u. When temperature is low enough (fc^T <C 
r, fiuo)i there is a well-defined boundary for the occupied 
states and empty states in the exploring terminal. With 
a change of the terminal voltage V3, the change of inci- 
dent electrons concentrate at a very small energy region, 
thus the differential conductance gives an excellent map- 
ping of the spectral function. On the other hand, this 
well-defined boundary for the occupied and empty states 
is destroyed when the temperature ksT ~ T, so that 
the conductance G 3 and the spectral function A(uj) have 
a large discrepancy. Naturally, if the third lead can be 
individually fixed at low temperature, then the spectral 



function can still be obtained from the differential con- 
ductance of the third lead regardless of the temperature 
in other parts of the system. In fact, the situation with 
low temperature for the third lead is always comparable 
with the above discussed low temperature regime. Let us 
discussed the realizability of the low temperature condi- 
tion fcfiT r, fkoo in the present technology. In an ex- 
periment, the characteristic frequency of phonon is about 
from 5meV to 50meV^^ and the coupling strength T is 
usually in the order of 100/ieU. But the temperature can 
reach 50mK in the present technology. So it should be 
achievable for the condition fc^T <C T, hujQ. In addition, 
in experiments the third weakly-coupled lead can be a 
probing STM tip, and then the coupling strength can be 
easily controlled by adjusting the distance between the 
STM tip and the SMT. 

At end, we make one more comment. Since we have 
used the same approximation as in Ref.(— ), the spectral 
function A{oj) is completely same as in their work. If 
to take a different approximation [e.g. as in Ref.(— )], 
the spectral function A(w), in particular the shape of the 
phonon-assisted side peaks, perhaps may vary somewhat. 
However, the perfect matching for the curves G3-V3 and 
A(uj)-lo still maintain as long as at low temperature and 
weak coupling conditions are met. 



V. CONCLUSIONS 

In summary, we study the transport behaviors of the 
three terminal SMT device coupled to a phonon mode. It 
is found that the intriguing characteristic of the phonon- 
assisted side peaks in the spectral function versus the 
energy can be directly observed from the differential con- 
ductance versus the voltage of the third weakly-coupled 
lead. In particular, not only the positions but also the 
shapes and the detailed structure of the side peaks of 
the spectral function can be perfectly mapped into the 
conductance if certain experimental conditions are met. 
Moreover we determine the conditions for this perfect 
mapping. The results exhibit that this mapping is ex- 
cellent at low temperature regardless of the bias and the 
level of the SMT. The mapping is destroyed at high tem- 
perature. 
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FIG. 2: (Color online) The dimensionless spectral function 
A(ui)/Ao vs. the energy ui at Ts = and the dimension- 
less differential conductance G3/G0 of the third lead vs. the 
voltage V3 at F3 = 0.01 for different so- Other parameters 
are taken as: ul = ur = 0, ujo = 5, T = 0.02, and A = 3. 
The unites Aq and Go are equal to 2/r and (2e 2 /h){T 3 /T) 
respectively. Notice that the three dotted curves almost over- 
lap perfectly with the three solid curves so that they almost 
cannot be seen in the figure. 



FIG. 1: Schematic diagram of a SMT coupled to three leads 
with tunnelling strength T a and the electron in the SMT is 
also coupled to a single-phonon mode. A gate electrode is 
capacitively attached to the SMT to tune the energy level of 
SMT. 
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FIG. 3: (Color online) The dimensionless spectral function 
A(ui)/Aa vs. the energy ui at r$ = and the dimensionless 
differential conductance G3/G0 of the third lead vs. the volt- 
age V3 at Ts = 0.01 for different £0, with ul = —ur = 3. 
The other parameters are same with Fig. 2. Notice that the 
three dotted curves almost overlap perfectly with the three 
solid curves so that they almost cannot be seen in the figure. 
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FIG. 4: (Color online) The dimensionless spectral function 
A(ui)/Ao vs. the energy u> and the dimensionless differential 
conductance G3 of the third lead vs. the voltage V3 for differ- 
ent coupling strengths F3. The parameters are taken as: uiq = 
5, T = 0.02, and A = 3. (a,c) and (b,d) are corresponding to 
ul = ur = and £0 = 0, and ul = — ur = 3 and eo = 6, 
respectively. The unites A = 2/T, Go = (2e 2 /^-)( r 3/ r ). and 
G' is determined by the equation / = J 



FIG. 5: (Color online) The dimensionless spectral function 
A(ui)/Aa vs. the energy ui at r$ = and the dimensionless 
differential conductance G3/G0 of the third lead vs. the volt- 
age V3 at Fa = 0.01 for different temperatures. The other 
parameters are ul = Ur = 0, ujq = 5, eo = 0, and A = 3. The 
unites A = 2/T and Go = (2e 2 /ft)(r 3 /r). 



